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ABSTRACT 

Linear stability theory is used to study the effect of crossflow on Gortler instability in 
incompressible boundary layers. The results cover a wide range of sweep angle, pressure gradient, 
and wall curvature parameters. It is shown that the crossflow stabilizes Gortler disturbances by 
reducing the maximum growth rate and shrinking the unstable band of spanwise wave numbers. On 
the other hand, the effect of concave wall curvature on crossflow instability is destabilizing. 
Calculations show that the changeover from Gortler to crossflow instabilities is a function of Gortler 
number, pressure gradient and sweep angle. The results demonstrate that Gortler instability may 
still be relevant in the transition process on swept wings even at large angles of sweep if the 
pressure gradient is sufficiently small. The influence of pressure gradient and sweep can be 
combined by defining a crossflow Reynolds number. Thus, the changeover from Gortler to crossflow 
instability takes place at some critical crossflow Reynolds number whose value increases with 
Gortler number. 
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1. Introduction 


A problem of great practical interest in aerospace applications is the drag reduction on aircraft 
wings, nacelles, etc. This has stimulated a number of research efforts on laminar-turbulent 
transition phenomenon for effective application of laminar flow control (LFC) which offers significant 
potential for reducing drag on aerodynamic surfaces. The supercritical wing is an example of 
advanced designs of laminar-flow airfoils. This wing has two concave regions on its lower surface 
and, therefore, when it is swept, the flow is subject to not only crossflow instability but also possibly 
the Gortler instability. The crossflow instability is observed on swept wings and leads to streamwise 
corotating vortices while Gortler instability is identified by the formation of pairs of counterrotating 
vortices with their axes aligned with the flow direction. 

In this paper, the effect of crossflow on Gortler instability is investigated. In the case of swept 
supercritical wing which supports both Gortler and crossflow instabilities a question arises as to 
which of these two vortex instabilities is significant for a given set of conditions. This question is 
also relevant for hypersonic flow ahead of engine inlet where concave curvature and some three- 
dimensionality are present either by design or inadvertently. Previous works addressing this 
problem are few and limited in scope and include theoretical works of Hall (1985), Collier and Malik 
(1987), Bassom and Hall (1991), and Blackaby and Choudhari (1993) and the experimental work of 
Kohama (1987) and more recently by Bippes (1994). 

Hall (1985) studied the effect of crossflow on Gortler vortices for flow over an infinitely long 
swept cylinder using asymptotic analysis. He concluded that the crossflow has a stabilizing effect on 
Gortler vortices and that for large angles of sweep the Gortler instability does not exist. For the 
concave region of a supercritical wing, Collier and Malik (1987) used Orr-Sommerfeld approach and 
showed that at low sweep angles the instability is of Gortler type and as the sweep angle increases 
the instability becomes of crossflow type. They also showed that concave surface curvature has a 
destabilizing influence on crossflow instability. Using hot-wire anemometry and smoke-wire flow 
visualization techniques, Kohama (1987) observed Gortler instability in the concave region in the 
absence of sweep and crossflow instability for 47° of sweep. The relation between Gortler and 
crossflow vortices was further studied by Bassom and Hall (1991). They found that Gortler vortices 
cannot exist on swept wings when the angle of sweep exceeds 20°. Blackaby and Choudhari (1993) 
studied inviscid instability of three-dimensional boundary layers over both concave and convex 
surfaces. They also investigated the relation of this problem to the stability of stratified shear flows. 
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Most of the aforementioned theoretical studies were based on asymptotic analyses in the large 
Gortler number limit and within the framework of inviscid linear theory. Moreover, they were 
restricted to a single value of the pressure-gradient parameter (fih = 0.5). In this paper, we present 
detailed results on the effect of crossflow on Gortler vortex instability by using viscous analysis. The 
governing mean flow equations and linear stability analysis are presented in the next section 
followed by the results in section 3. 

2. Problem Formulation 

We consider a three-dimensional boundary-layer flow over a concave surface whose constant 
radius of curvature r* = l/x t . The stream wise, wall-normal and the spanwise coordinates are 
denoted as x = Xf( 0 , y = YU 0 and z = Z/t 0 , respectively (y = 0 denotes the wall), where the length 
scale t Q will be prescribed later. Let the x,y and z components of the velocity and pressure be given 
by 


[u\v\w 1 ) = U e {U(x,y) + u(x,y,z,t),V(x,y)+v(x,y,z,t),W + w{x,y,z,t)) (1) 

p f = pU' (P + p(x, y, z, t)) (2) 

where the superscript t represents a dimensional quantity and U e is the velocity scale. Here U, V 
and W are mean-flow components whereas u, v, w represent the perturbation velocity components in 
x, y, z directions, respectively. Similarly, P and p represent the mean and perturbation pressures. 

The mean flow investigated in this study is the Falkner-Skan-Cooke (FSC) boundary layer. We 
consider the chordwise potential flow velocity distribution: 

U„ = cX m 

where the exponent m is related to the Hartree pressure gradient parameter p h = 2m/(m + 1) which 
defines the wedge half angle (n / 2)p h . The spanwise potential flow velocity component, W_, parallel 
to the leading edge is constant. With the similarity variable J] defined by rj = YU where 

u w 

( = (vX/UJ U2 , vthe kinematic viscosity, U = -=y-f’, W = — r-g, the self-similar three-dimensional 

u e LJ € 

boundary layer equations are given by: 

f" + ff" + m(l- /"' 2 j = 0 (3) 

2 
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( 4 ) 


where the primes denote differentiation with respect to rj. 

We assume that the Reynolds number is large and that the radius of curvature is much larger 
than the boundary-layer thickness, S (i.e., k^S « 1). In this case, if k = kU 0 , the equations 
governing the perturbation quantities are 


du TT du BU TT du dU da , TT „ » dp 1 _2 « 

— +U— + u — + V— + v— + W— + k(Uv + Vu) + V u = 0 
dt dx dx dy dy dz dx R 


(5) 


du TT du dV Tr dv dV TI . du n dp l n2 „ 

— + U— + U— + V— + V— + W- — 2kUu+-£--—V v = 0 
dt dx dx dy dy dz dy R 


( 6 ) 


dw TT dw — ~ dw „ T dw dp 1 _ 2 n 

— + U— + V— + W— +-zr~— Vw = 0 
dt dx dy dz dz R 


(7) 


where 

The boundary conditions are 

u-v-w-0 at y = 0 and u -> 0, v -> 0, w -» 0, as y -> <» (9) 

Here, t 0 and U e are constants so as to define Reynolds number R as 

_ UJo 
v 

where the length scale = ^jvX 0 /U e , X 0 being the location (dimensional) of a reference streamwise 
station. In flows where Gortler vortex phenomenon constitutes the dominant primary instability, 
Gortler number is generally of 0(20) where G is defined as 

G = R^e (10) 


du do dw - 


(8) 


v 2 = 


dx 2 + dy 2 + dz 2 +K dy' 
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In the present paper, we solve the viscous linear stability equations under the parallel-flow 
approximation. While such an approximation is questionable for Gortler vortices at 0(1) Gortler and 
wave numbers, it is less so for high Gortler numbers. The quasi-parallel approach, despite its limi- 
tations, provides a useful tool to study the link between Gortler and crossflow disturbances, includ- 
ing the effect of viscosity. Parabolized stability equations (Bertolotti et al. 1992, Malik and Li 1993) 
represent a more appropriate model for this problem, but the local approach adopted here allows 
parametric studies to be performed more efficiently. Except for flat walls ( G = 0) where only 
crossflow instability is relevant and for which quasi-parallel assumption has been commonly used, 
we will consider large Gdrtler numbers (G £ 15) so that the error introduced by the above 
approximation is minimized. For large Gdrtler numbers, the present results are at least as valid as 
the previous asymptotic results with the added effect of viscosity. 

Under the parallel flow approximation (U = U(y), V = 0,W = W(y)) allows stationary 
disturbances of the form 


<t>(x,y,z ) = 0(y)e' (<tt+w +c.c. 

(id 

where <j> = [u,v,w,p] T and c.c. stands for complex conjugate. Substituting Eq. 
yields the following ordinary differential equations 

(11) into Eqs. (5H8) 

ju + iap = 0 

(12) 

-2tdJu + $3 + — = 0 

dy 

(13) 

%w + = 0 

(14) 

iau + ^- + ipw + kv = 0 

dy 

(15) 

where 


S = iaU+ipW--^ - -a 2 -p 2 

R[ dy 2 dy 


The boundary conditions for the above equations are 



u(0) = 5(0) = i£(0) = 0 


(16) 
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« -» 0, 0 -> 0, w -> 0 as y-> » 


(17) 


Equations (12H15) along with homogeneous boundary conditions (16H17) constitute an eigenvalue 
problem 


a = a(p) (18) 

which is solved by fourth-order accurate compact difference scheme described by Malik et al. (1982). 
As noted in Malik and Li (1993) and Masad and Malik (1994), the dominant curvature term for 
three-dimensional boundary-layer stability is the 2 kUu term in the normal momentum equation. An 
order of magnitude analysis of the Gortler vortex problem in two-dimensional flows leads to the same 
conclusion (Hall 1983). 

3. Results 

The results presented in this paper are obtained by solving the linear stability equations (12-17) 
for a Reynolds number of R = U e £ 0 /v= 1000. While the Gflrtler instability depends only on the 
value of the Gortler number (Eq. (10)), crossflow instability depends also on the value of Reynolds 
number, but at this high value of Reynolds number the dependence is weak since crossflow 
instability is basically inviscid in nature. Figure 1 (for = 0.5) shows the effect of sweep angle A in 
radians A = tan’^W^/C/^), on spatial growth rate (scaled with / 0 ) of stationary disturbances for 
different Gortler numbers. The corresponding streamwise wave numbers are shown in Fig. 2. In the 
absence of wall curvature (Fig. la) where the instability is of crossflow type the growth rate 
increases with the increase in the sweep angle. As Gfirtler number increases (Figs, lb, c, d) two 
distinct features of the growth rate are observed; at low sweep angles the instability is of Gortler 
type and as the sweep angle increases beyond a certain limit the instability becomes of crossflow 
type. This sweep angle limit corresponding to the changeover from Gortler to crossflow instabilities 
is dependent on Gortler number and as will be seen below, on the pressure-gradient parameter . 
As shown in Fig. 1 the "changeover" sweep angle increases with the increase in Gortler number. For 
G - 50, the changeover occurs at approximately A = 0.4 radians (23°). The corresponding value 
reported by Bassom and Hall (1991) for the large Gortler number inviscid limit is A = 20°. Thus our 
viscous calculation using j3 h = 0.5, G = 50 yields a result close to that obtained by the inviscid theory. 

In the sweep angle range over which Gortler instability is present, we notice (see Figs. 1 b,c,d) 
that the maximum growth rate decreases with the increase in the sweep angle. Moreover, the 
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unstable band of spanwise wave numbers shrinks with the increase in crossflow; thus, the effect of 
crossflow on Gortler disturbance is a stabilizing one. It is also to be noted that in the inviscid Gortler 
problem, only the left branch of the growth curve is present. The right branch which gives rise to a 
neutral point at high wave numbers is due to the effect of viscosity. This right neutral point moves 
to lower wave numbers when sweep angle increases as noted above. 

Another observation is made here with regards to the influence of wall curvature on crossflow 
instability. The results shown in Fig. 1 are replotted in Fig. 3 for six different sweep angles. It is 
seen that concave wall curvature enhances crossflow instability (see Figs. 3 c,d,e,f), as was earlier 
shown by Collier and Malik (1987). In contrast, convex wall curvature stabilizes crossflow 
disturbances (Malik and Balakumar 1993, Masad and Malik 1994). These effects are also consistent 
with the inviscid results of Blackaby and Choudhari (1993). 

An important parameter relevant in three-dimensional boundary-layer flows is the crossflow 
Reynolds number defined as 


( 19 ) 

where U e is the magnitude of the maximum crossflow velocity (in a direction which is at right angle 
to the local inviscid streamline) within the boundary layer and S 01 is the boundary-layer thickness 
within which the crossflow velocity drops to 10 percent of U c . 

In Fig. 4, the variation of crossflow Reynolds number with sweep angle is plotted for different 
P h . It is clear that the crossflow Reynolds number increases pressure gradient for any sweep angle. 
Therefore, the crossflow instability is expected to be enhanced with pressure gradient. The results in 
Figs. 1-3 were obtained for a single value of the pressure gradient ( p h = 0.5). The results for 
different values of pressure gradient parameter are presented in Fig. 5. The effect of pressure 
gradient is seen to influence Gdrtler and crossflow instabilities in two different ways. As the 
pressure gradient becomes more favorable, the maximum Gortler vortex growth rate decreases and 
the unstable band of spanwise wave numbers shrinks. In contrast, the crossflow instability 
increases with the pressure gradient. Thus, the increase in p h destabilizes crossflow vortices and 
stabilizes Gortler vortex instability. This is a direct result of increase in the crossflow velocity (and 
consequently R cf ) with pressure gradient. 

For any given sweep angle, results such as those in Fig. 5 enable us to identify the pressure 
gradient at which Gortler instability is destroyed in favor of crossflow instability. Alternatively, for 
any given pressure gradient the sweep angle corresponding to the changeover from Gortler to 
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crossflow instabilities can be determined. When the maximum growth rate is calculated for different 
values of pressure gradient and sweep angles and the results are plotted in the manner shown in 
Fig. 6, the following observations can be made: 

(1) Each curve corresponding to a given sweep angle has, when applicable, two branches of opposite 
slopes. The lower branch (negative slope) corresponds to Gortler instability while the upper 
branch corresponds to crossflow instability. 

(2) For small sweep angles (A *= 0.1 radians) the instability is of Gortler type for any pressure 
gradient parameter P h . 

(3) For large sweep angles Gortler instability is restricted to small pressure gradients. 

Thus, a reduction in the pressure gradient renders instability of the Gortler type. For example, 
for G = 15 (Fig. 6a) Gortler instability is the dominant mechanism for pressure gradients P h < 0.25 
and sweep angles greater than 0.3 whereas for G = 30 (Fig. 6b) this occurs at higher pressure 
gradients. Thus, Gortler instability may be the cause of transition in incompressible flow even at 
high angles of sweep if the pressure gradient is sufficiently small but nonzero. For zero-pressure 
gradient, R cf = 0 for any sweep angle and there is no crossflow instability. 

In three-dimensional boundary-layer transition, crossflow Reynolds number, R c f, is used as a 
governing parameter in correlating the transition onset. In this study we look into a possible 
correlation for the changeover from Gortler to crossflow instability in terms of R c f. Figure 7 shows 
the variation of maximum growth rate with R cf for different pressure gradient levels while Fig. 8 
shows the results for fixed sweep angles. The R cf double values (e.g.. Fig. 7b for p h = 0.5) are a 
result of the R cf variation with sweep angle shown in Fig. 4 as A increases beyond a value of 0.76 
radians. In Figs. 7a, b for any fixed pressure gradient, the minimum value of the maximum growth 
rate, when applicable, corresponds to the changeover from Gortler (left branch) to crossflow (right 
branch) instability. The corresponding R cf denoted as R cfc (the changeover R cf ) is seen to be a 
function of pressure gradient, sweep and the Gortler number. However, the strong dependence of 
R c f c is on Gortler number. A comparison of Figs. 7(a) and 8(a) for G = 15 shows that R cfc ~ 70-100 
while a similar comparison of Figs. 7(b) and 8(b) for G =30 shows that the changeover takes place at 
R c f c ~ 150-200. Thus R c f is an important parameter for determining the changeover from Gortler to 
crossflow instability but its value increases with Gortler number. 
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The above results represent an important guide for the design of advanced LFC wings. 
However, it is cautioned that they are based upon local theory. Clearly, whether crossflow can 
destroy Gortler vortices depends upon nonlinear effects and the strength of the vortex which is 
related to the upstream history (distribution of wall curvature and pressure gradient for a given 
leading-edge sweep) and this evolution is best described by nonlinear parabolized stability equations. 


4. Conclusions 

An investigation of the relation between Gortler and crossflow instabilities of incompressible, 
three-dimensional boundary-layer flow is carried out using linear stability theory. The results cover 
a wide range of sweep angles, pressure gradients and wall curvature. It is shown that Gortler 
instability may persist even at large sweep angles when the pressure gradient is small. The results 
demonstrate the stabilizing effect of crossflow on Gortler disturbances and the destabilizing effect of 
concave wall curvature on crossflow disturbances. It is shown that the changeover from Gortler to 
crossflow instabilities is governed by crossflow Reynolds number. The value of crossflow Reynolds 
number at which the instabilities are switched increases with increasing Gortler number. 
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Figure Captions 


Figure 1. Effect of sweep angle A (radians) on spatial growth rate of stationary disturbances in 
Falkner-Skan-Cooke (FSC) boundary layer for different values of Gortler numbers. 

Figure 2. Variation of streamwise wave number with spanwise wave number for different sweep 
angles and Gortler numbers. 

Figure 3. Effect of wall curvature on spatial growth rate of Gortler and crossflow disturbances in 
FSC boundary layer for different values of sweep angle A. 

Figure 4. Variation of crossflow Reynolds number with pressure gradient and sweep. 

Figure 5. Effect of pressure gradient on spatial growth rate of stationary crossflow and Gortler 
disturbances of FSC boundary layers. 

Figure 6. Effect of sweep angle and pressure gradient on maximum growth rate of stationary 
crossflow and Gortler disturbances of FSC boundary layer for two different Gortler numbers. 

Figure 7. Variation of maximum growth rate with R cf for fixed P h . 

Figure 8. Variation of maximum growth rate with R c f for fixed A. 
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Figure 2. Variation of streamwise wave number with spanwise wave number for different sweep 
angles and Gortler numbers. 
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Figure 6. Effect of sweep angle and pressure gradient on maximum growth rate of stationary 
crossflow and Gortler disturbances of FSC boundary layer for two different Gortler numbers. 
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